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The liquid-gas phase transition in hot asymmetric nuclear matter is investigated within relativis- 
tic mean-field model using the density dependence of nuclear symmetry energy constrained firom 
the measured neutron skin thickness of finite nuclei. We find symmetry energy has a significant 
influence on several features of liquid-gas phase transition. The boundary and area of the liquid-gas 
coexistence region, the maximal isospin asymmetry and the critical values of pressure and isospin 
asymmetry all of which systematically increase with increasing softness in the density dependence of 
symmetry energy. The critical temperature below which the liquid-gas mixed phase exists is found 
higher for a softer symmetry energy. 
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The possible occurrence of liquid-gas phase (LGP) 
transition in intermediate energy heavy ion collisions us- 
ing neutron-rich stable and future radioactive beams pro- 
vides a rather unique tool to probe hot and dense phases 
of highly asymmetric nuclear matter. Collisions exper- 
iments [l|, 0] with stable heavy nuclei at intermediate 
energy do indicate theoretically predicted Q features of 
liquid-gas phase transition where the hot and compressed 
nucleus produced expands and fragments into several in- 
termediate mass fragments (high-density liquid phase) 
and light particles and nucleons (low-density gas phase). 

The early theoretical studies of the thermodynamic 
properties of liquid-gas phase transition ^-(Tj] are mostly 
confined to symmetric nuclear matter that employed the 
quite well predicted [8l-[l0| behavior of the symmetric nu- 
clear matter equation of state (EOS). One of the major 
ingredients in studies of asymmetric nuclear matter re- 
quire knowledge of the density dependence of symme- 
try energy Egy^nip) [TTI - [T3| . Unfortunately, the model 
predictions of -Bsyin(p) even for nuclear matter at zero 
temperature are extremely diverse [14i] . Only at the nu- 
clear saturation density po w 0.16 fm"-^ the value of 
E{po, T = 0) = 32 ± 4 MeV has been well constrained. 

Recently some progress has been achieved by consis- 
tently constraining the symmetry energy of cold neutron- 
rich matter near normal matter density from analysis of 
isospin diffusion [Tl| - [T3| and isoscaling [l^ data in inter- 
mediate energy heavy ion collisions and from the study 
of neutron skin thickness of several nuclei [l^ [13] ■ While 
knowledge of symmetry energy £^sym(p, T') at finite tem- 
perature in particular has received little attention [l8l - [20j 
that is crucial for a proper understanding of the features 
of LGP transition in hot asymmetric nuclear matter. 
In fact, new qualitative features are expected when an 
asymmetric nuclear system with two conserved charges, 
baryon number and third component of isospin, under- 
goes a LGP change which has been suggested to be of 
second order [2l|. Most previous studies of LGP tran- 
sition [2l| - |23j relied on model predictions of symmetry 
energy Esyin{p,T) with no or minimal contact with the 



available experimental data. Thus to understand better 
the features of LGP transition in hot asymmetric nu- 
clear matter, it is imperative to employ the asymmetric 
nuclear EOS that has been constrained from analysis of 
skin thickness data of several nuclei or from isospin 
diffusion/scaling data [H, [1^. Such an investigation is 
particularly useful as future experiments with radioac- 
tive ion beams with large neutron-proton asymmetries 
can be used to explore 0, HE] symmetry energy effects 
on liquid-gas phase transition. 

In this letter, we study the effects of constrained sym- 
metry energy [l7[ on the thermodynamic properties of 
LGP in hot neutron-rich nuclear matter within relativis- 
tic mean field (RMF) models [26l|. For this purpose 
we use two accurately calibrated models: NL3 [23] and 
FSUGold [il], that was obtained by fitting the model 
parameters to certain ground state properties of finite 
nuclei. The interaction Lagrangian density in the non- 
linear RMF model is given by 
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which includes a isospin doublet nucleon field (ip) in- 
teracting via exchange of isoscalar-scalar sigma ((/)), 
isoscalar-vector omega (V^), isovector-vector rho (b'^) 
meson fields and the photon (v4^) field. The nonlinear 
sigma meson couplings (k, A) soften the symmetric nu- 
clear matter EOS at around po, while its high density 
part is softened by the self-interactions (C) for the omega 
meson field. 

For the original NL3 set with ^ = A^, =0, the sat- 
uration of symmetric nuclear matter occurs at a Fermi 
momentum of kp — 1.30 fm^^ with a binding energy 
B/A Ri 16.3 MeV and an incompressibility of Kq = 271 
MeV. The original FSUGold with two additional 
couphngs C = 0.06 and A„ = 0.03, with Kq = 230 MeV 
produces a soft symmetric and asymmetric nuclear mat- 
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ter EOS. To study the effect of symmetric nuclear EOS 
(eg. incompressibility Kq) on tire symmetry energy, we 
fiave extended p/fl] tlie original NL3 Lagrangian to in- 
clude the isovector coupling A^, which is then varied along 
with gp in both NL3 and FSUGold to generate various 
Esymip)- AH combinations of At, and gp are adjusted 
to a constant Esy^{p,T = 0) = 25.67 (26.00) for the 
NL3 (FSUGold) at an average density p corresponding 
to kp = 1.15 fm"'^ where the binding energy of ^"^^Pb 
is reproduced. Thus the additional couplings provides 
an efficient way to change in a controlled manner the 
density dependence of nuclear symmetry energy without 
compromising the success of the model. 

The model parameter (Au, gp) set is then varied to ex- 
plore i?syin (p) effects on the liquid-gas phase transition in 
hot asymmetric nuclear matter. For the present study we 
use A,; = 0.0 — 0.03 since the resulting symmetry energies 
and their slopes and curvatures are in reasonable agree- 
ment with that extracted from neutron skin thickness of 
several nuclei as well as the isoscaling and isospin diffu- 
sion data [13] • It may be also noted that with increasing 
A„ the density dependence of symmetry ener gy b ecomes 
softer in both the NL3 and FSUGold models [ill. While 
at a finite A,; the symmetry energy Esym{p,T = 0) is 
found to be particularly stiff in FSUGold than in the 
NL3 parameter sets at densities p ^ 1.5/3o- 

At finite temperature and density the energy density 
£ can be readily obtained from the thermodynamical po- 
tential Vt 211 as 
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where E*{k) ~ ^/k"^ + ra*"^ is the effective energy. The 
distribution function for nucleon and antinucleon (re- 
ferred to as ± sign) 
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where the effective chemical potential for neutron and 
proton is expressed a.suq = pq—gyVozLgpbo/2. The chem- 
ical potentials can be determined from the conserved 
baryon and isospin densities: 
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FIG. 1: Density dependence of nuclear symmetry energy at 
temperatures T = 0, 5, 10, 15 MeV in the NL3 (left panel) and 
FSUGold set (right panel) with couplings A,, — 0.0 and 0.03. 



hot neutron-rich nuclear matter can be expressed in the 
parabolic form: 



E{p,T,a) ^ E{p,T,a ^ 0)+E,y^{p,Ty 
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where the neutron-proton asymmetry is a = (p„ — pp)/p. 
The density and temperature dependence of symmetry 
energy can be estimated from Esymip^T) w E{p,T,a = 
1) - E{p,T,a = 0). This implies that E^ymip^T) is the 
energy required to convert all the protons in symmet- 
ric matter to neutrons. Figure [T] shows the density de- 
pendence of nuclear symmetry energy at temperatures 
T = 0, 5, 10, 15 MeV in the NL3 (left panel) and FSUGold 
(right panel) sets. For all choices of A,, the symmetry en- 
ergy decreases with increasing temperature especially at 
small densities p po that is entirely due to the decrease 
in the kinetic energy contribution. For A^, = 0.0 (0.03) 
the density dependence of £^sym (p, T) at all temperatures 
exhibits a systematic trend of small (large) value at sub- 
saturation densities and a large (small) value at supranor- 
mal densities resulting in an overall stiffer (softer) asym- 
metric nuclear matter EOS. 

The above described models can now be used to study 
LGP in hot asymmetric nuclear matter. The system is 
stable against LGP separation if its free energy F is lower 
than the coexisting liquid {L) and gas (G) phases, i.e. 
F{T,p) < {1- \)F^{T,p^) + \F^{T,p'^) withp= (1- 
X)p^ + Xp'^ where < A < 1 and A = /V being the 
fraction of the total volume occupied by the gas_phase. 
The stability condition implies the inequalities [21i] : 
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where Gq{k) = K(fc)]^ - [ng(fc)]_. 

As in the zero temperature case, several model stud- 
ies [H [1^, m, [13 have indicated that the EOS for 
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FIG. 2: Pressure as a function of density at temperature T = 
10 MeV for various isospin asymmetry a in the original NL3 
set pTIl with A„ — 0.0 (top panel) and the original FSUGold 
set [2^ with A„ = 0.03 (right panel). The dotted curves refer 
to unstable single phase while the solid curves refer to stable 
matter; see text for details. 



The first inequality indicates mechanical stability which 
means a system at positive isothermal compressibility re- 
mains stable at all densities. The second inequality stems 
from chemical instability which shows that energy is re- 
quired to change the concentration in a stable system 
while maintaining temperature and pressure fixed. If one 
of these conditions get violated, a system with two phases 
is energetically favorable. The two phase coexistence is 
governed by the Gibbs's criteria for equal chemical po- 
tentials and pressures in the two phases with different 
densities but at the same temperature: 

M^(T,/) = M^(T,p«) (g = n,p), (9) 
F^(T,/) = F^(r,pG). (10) 

Figure [5] shows the pressure as a function of nucleon 
density at a fixed temperature T = 10 MeV with different 
values of asymmetry a in the original NL3 and FSUGold 
sets. Below a critical value of asymmetry a, the pressure 
is seen (dotted curves) to decrease with density result- 
ing in negative incompressibility and thereby a mechani- 
cally unstable system. The stable two-phase (liquid-gas) 
configuration at each density is obtained from Maxwell 
construction (solid lines). Analogues to intermediate en- 
ergy heavy- ion collisions [l|, 0| when the hot matter in the 
high density (liquid) phase expands it enters the coexis- 
tence LGP where the pressure decreases at a fixed a 7^ 
for the two-component asymmetric matter. Whereas, for 
symmetric nuclear matter at a = the pressure remains 
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FIG. 3: Chemical potential isobars as a function of isospin 
asymmetry a at temperature T = 10 MeV for NL3 (left panel) 
and FSUGold (right panel) with different A„ couphngs. The 
geometrical construction used to obtain the isospin asymme- 
tries and chemical potentials in the two coexisting phases is 
also shown. 



constant at all densities. Finally the system leaves the co- 
existence region and vaporizes into the low density (gas) 
phase. Of particular interest here is the symmetry en- 
ergy effects on the isotherms. It is clearly seen that in 
contrast to the original NL3 with A„ = 0, the softer 
£^sym(p) in the original FSUGold with A^, = 0.03 [13 en- 
forces the onset of pure liquid phase to a higher density 
resulting in a wider coexistence region for each asymme- 
try a. Moreover, the critical pressure Pc above which the 
mixed liquid-gas phase vanishes is seen larger for this soft 
FSUGold set; a detailed discussion of which is presented 
below. 

The details of chemical evolution for the LGP transi- 
tion is depicted in Fig. [3] where the neutron and proton 
chemical potentials are shown as a function of isospin 
asymmetry a at a fixed T ~ 10 MeV and pressure 
P = 0.11 MeV/fm^ for the NL3 (left panels) and FSUG- 
old (right panels) at various A„ values. As usual, the 
bare nucleon mass has been subtracted from the chemi- 
cal potentials. At fixed pressure and A^,, the solutions of 
the Gibbs conditions ^ and ([TU|) for phase equilibrium 
form the edges of a rectangle and can be found by geo- 
metrical construction as shown in Fig. [3l At each A^,, the 
two different values of a defines the high density liquid 
phase boundary (with small a = ai(r, P)) and the low 
density gas phase boundary (with large a = a2(P, P)). 
From the figure it is evident that the symmetry energy 
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dependence of A„ in NL3 and FSUGold [17|, l25| leads 
to different phase boundaries ai{T,P) and a2{T,P) and 
hence should predict different thermodynamic properties 
for the LGP transition. 

As the pressure increases the system encounters a 
critical pressure P^. beyond which the matter is stable 
but below which the second inequality ^ gets violated 
and the system becomes chemically unstable. The crit- 
ical pressure Pc is determined by the inflection point 
{dfi/da)T,p^ = {d'^fi/da^)T,p^ = 0. The disappearance 
of chemical instability at Pc results in the neutron (pro- 
ton) chemical potential to decrease (increase) with de- 
creasing asymmetry a. Figure [3] also shows the chemical 
potential isobars at the critical pressure (dashed lines). 
The rectangle from Gibbs condition then collapses into a 
line vertical at a = ac- Correspondingly, {P^ ac) defines 
the critical point at a given temperature that refers to 
the upper boundary of instability with respect to pres- 
sure variation. Note at T = 10 MeV, the critical val- 
ues {Pc,ac) at Ay = 0.0,0.02,0.03 are respectively at 
(0.210, 0.652), (0.276, 0.741), (0.331, 0.797) for the NL3 
set and at (0.209, 0.638), (0.266, 0.725), (0.303, 0.789) 
for the FSUGold set. Interestingly, we also find at a fi- 
nite temperature the stiffness of symmetry energy has a 
significant infiuence on the phase-separation boundaries 
of LGP transition In general, a softer symmetry en- 
ergy (larger A^,) gives systematically larger critical pres- 
sure and an enhanced asymmetry in the system. More- 
over at a finite a, the relatively softer symmetry energy 
Esyrn{p,T = 0) in the NL3 compared to FSUGold [iTl 
translates to a larger critical pressures and asymmetry 
for the LGP transition. 

All the pairs of solutions of Gibbs conditions, ai{T, P) 
and a2 (T, P) , form the phase-separation boundary or the 
binodal surface. In Fig. |4] we show the section of the 
binodal surface under isothermal compression of asym- 
metric nuclear matter at T = 10 MeV in the NL3 (top 
panel) and FSUGold (bottom panel). As expected the 
point of equal concentration (EC) corresponding to sym- 
metric nuclear matter is independent of A„. The criti- 
cal point (CP) and EC divide the binodal section into 
two branches. One branch is the high-density (liquid) 
phase that is less asymmetric while the other branch 
corresponds to the more asymmetric low-density (gas) 
phase. Thus the matter on the left (right) of the binodal 
surface represents stable liquid (gas) phase. It is clearly 
seen here that the critical point {Pc,ac) depends on the 
density dependence of the symmetry energy associated 
with different A„ values. 

We also indicate on the binodal surface the maximal 
isospin asymmetry (MA), auA, of the system. Thus 
more neutron-rich matter on the right side of the sur- 
face when compressed/expanded at fixed a will never 
encounter a coexistence phase. Note here the maximal 
asymmetry is also quite sensitive to A^, i.e. on i?sym(p, T). 
Such effects found in the present study should have strong 
infiuence on the experimentally observed isospin distilla- 
tion phenomena |31j where the gas phase is more neutron- 
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FIG. 4: The section of binodal surface at temperature T = 10 
MeV in NL3 (top panel) and FSUGold (bottom panel) with 
different At, couphngs. The critical point (CP), the point of 
equal concentration (EC), and the maximal asymmetry (MA) 
are indicated. 



rich (large n/p ratio) compared to the more asymmetric 
liquid phase. However for pressures P > 0.10 MeV/fm^ 
the magnitude of isospin distillation is more sensitive to 
the symmetry energy used. 

A new feature for LGP transition in asymmetric sys- 
tem, refereed to as retrograde condensation [2l|, arises 
when a nucleon gas prepared at an asymmetry etc < 
a < ttMA is compressed at fixed total a. The matter re- 
mains mechanically stable but chemically unstable. Thus 
a coexisting liquid phase emerges which finally vanishes 
when the system leaves the binodal surface as a pure gas. 
As the extent Aa — aMA — etc is found to decrease for 
softer symmetry energy with higher A^,, the possibility 
of such unique-phase condensation phenomena also be- 
comes minimal. 

The present study clearly suggests that for liquid-gas 
phase transition in hot asymmetric nuclear matter, the 
critical values of pressure and isospin asymmetry, the 
maximal asymmetry and the area and shape of the bin- 
odal surface are quite sensitive to the density dependence 
of symmetry energy with a stiffer symmetry energy leads 
to consistently smaller values of these thermodynamic 
variables. 

The existence of critical isospin asymmetry parameter 
ac at a given temperature indicates that for a > ac the 
system will not change completely into the liquid phase. 
Conversely, this suggests that at a fixed a there exists a 
critical temperature Tc beyond which the system can only 
be in the gas phase at all pressures. In Fig. [S]we present 
Tc as a function of a in the FSUGold set for different 
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FIG. 5: The critical temperature Tc versus isospin asymmetry 
a for different A„ in the FSUGold set. 

couplings A„. For symmetric nuclear matter (a — 0), 
the critical temperature for LGP transition in this model 
is Tc = 14.7 MeV. With increasing asymmetry a ^ 0.6, 
Tc decreases rapidly. A softer density dependence in sym- 
metry energy (larger Au) shows the coexisting liquid-gas 



phase can prevail for larger values of Tc- We find that for 
the soft symmetry energy (A^, = 0.03) even pure neutron 
matter (a = 1) can exhibit LGP transition at T < Tc = 2 
MeV. While the stiffest symmetry energy [Ky — 0) at 
a > 0.9 predicts that the matter can only be in the pure 
gas phase at all temperatures. 

In summary the effects of isospin symmetry interaction 
on the liquid-gas phase transition in hot neutron-rich nu- 
clear matter is investigated. For this we have used the 
two accurately calibrated relativistic mean field models, 
the NL3 [131 and the FSUGold [2i| wherein the density 
dependence of nuclear symmetry energy at zero tempera- 
ture has been constrained within a limited range by neu- 
tron skin thickness data of several atomic nuclei. We 
find considerable sensitivity of the symmetry energy on 
the features of phase transition. Softer symmetry en- 
ergies give progressively larger phase-separation bound- 
aries with higher critical values for pressure and isospin 
asymmetry as well as maximal asymmetries. At a given 
asymmetry we find the critical temperature for the exis- 
tence of the mixed liquid-gas phase increases with softer 
symmetry energy and predicts the possible occurrence of 
even an unstable pure neutron matter at finite tempera- 
tures. 
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